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Abstract. Measurements of the cosmological density parameter (fi) 
using techniques that exploit the gravity-induced motions of galaxies 
yield, in linear perturbation theory, the degenerate parameter combina- 
tion (3 = where the linear bias parameter b is the ratio of the 
fluctuation amplitudes of the galaxy and mass distributions. However, 
the relation between the mass and the galaxy density fields depends on 
the complex physics of galaxy formation, and it can in general be non- 
linear, non-local, and stochastic. There is a growing consensus that the 
one-parameter model for bias is oversimplified. This leaves us with the 
following question: What is the quantity that is actually being measured 
by the different techniques? In order to address this question, we present 
estimates of /3 from galaxy distributions constructed by applying a variety 
of biased galaxy formation models to cosmological N-body simulations. 
We compare the values of /3 estimated using two different techniques: the 
anisotropy of the redshift-space power spectrum and an idealized version 
of the POTENT method. In most cases, we find that the bias factor 
b = Q^-^/ (3 derived from redshift-space anisotropy or POTENT is similar 
to the large-scale value of bfj{R) defined by the ratio of rms fluctuation 
amplitudes of the galaxy and mass distributions. However, non-linearity 
of bias and residual effects of non-linear gravitational evolution both in- 
fluence (5 estimates at the 10-20% level. 



1. Analysis 

We estimate (3 from a variety of simulated galaxy distributions, created by apply- 
ing simple local and deterministic biasing algorithms to the mass distributions 
obtained from cosmological N-body simulations. 

The cosmological models that we use are: (1) = 1.0, CDM model 
with a tilted power spectrum designed to simultaneously satisfy both COBE 
and cluster normalization constraints. (2) O = 0.4, CDM, cluster normalized 
model. (3) = 0.2, CDM, cluster normalized model. Cluster normalization 
requires that a^^nSl^-^ = 0.55 (White, Efstathiou & Prenk 1993). 

The biasing algorithms that we use are: (1) Semi-analytic: An empirical 
bias prescription derived by Narayanan et al. (1999, in prep.) which character- 
izes the relation between the galaxy and mass density fields in the semi-analytic 
galaxy formation models of Benson et al. (1999). (2) Sqrt-Exp. (Square-root 
Exponential): A biasing prescription in which {l + 5g) oc •\/(T+"5^e"(^"*'^™'*. (3) 
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Power-law: A biasing prescription in which (1 + Sg) oc (1 + Sm)"- Here, (1 + Sg) 
and (1 + Sm) are the galaxy and mass over-densities, respectively. 

All the bias models are designed to yield galaxy distributions with an rms 
fluctuation, in top-hat spheres of radius 12/i~^Mpc, of ai2 ~ 0.7. The values of 
n^-^/b„{12) for the galaxy distributions are 0.6, 0.58, and 0.56 for Q = 1.0, 0.4, 
and 0.2, respectively. All the results have been averaged over four independent 
realizations of the cosmological models. The simulation volumes are periodic 
cubes, with a box size of mh-^Mpc. 

The first method that we use to estimate P involves the anisotropy of the 
rcdshift-space power spectrum. Kaiser (1987) and Cole, Fisher Sz Weinberg 
(1994) showed that in linear theory, and under the plane-parallel approximation, 
the multipole moments of the redshift-space power spectrum may be used to 
estimate /3. Here we employ two different /? estimators: P^{k)/P^{k), the 
ratio of the angle-averaged redshift-space power spectrum (monopolc) to the 
real-space power spectrum, and P2{k)/PQ{k), the ratio of the quadrupole and 
monopole moments of the redshift-space power spectrum. These ratios are, in 

principle, measurable, and they are related to (3 as follows: -piq^ = 1 + i/3 + 

|/?^ ^ = (|/3 + f/?^)/(l + §/? + |/3^)- These ratios yield an estimate of 
f3 at each wavenumber k, and the scale dependence of this estimate is caused 
mainly by non-linearity in the peculiar velocity field. It is possible to obtain a 
global estimate of /? by modeling this non-linearity. We use two such models 
in our analysis: (1) Exponential velocity distribution model (Cole, Fisher & 
Weinberg 1995), where we assume that the galaxies have uncorrelated small 
scale peculiar velocities which are drawn from an exponential distribution. We 
use this model in our analysis of P^{k)/P^{k) (the results of which are shown in 
Fig. 1). (2) Empirical model (Hatton & Cole 1999), where the scale dependence 
of P2{k)/PQ{k) is found empirically by examining a large number of N-body 
simulations spanning a broad range of parameter space. 

The second method that we use to estimate (3 involves a direct comparison 
between the mass and galaxy density fields. The POTENT method (Bertschinger 
&; Dekel 1989) reconstructs the full three dimensional peculiar velocity field v 
from the observed radial peculiar velocity field, allowing a measurement of [3 
directly from the linear theory relation V • v = —[3H5g. We estimate (3 using an 
idealized POTENT method, in which we have perfectly reconstructed the three 
dimensional peculiar velocity field. Wc obtain the value of [3 by fitting a line of 
the form y = ax + c to the — (V • v) vs. 5g relation. We perform this fit in the 
region: —0.5 < 5g < 0.5. By repeating this procedure for different smoothing of 
the velocity and density fields, we can measure /3 as a function of scale (see Fig. 
!)• 

2. Conclusions 

Figure 1 compares (3 estimates from the two methods described previously to the 
function Pa{R) = ^^°'^/6o-- Here, ba{R) is the bias function defined as bcr{R) = 
(Tg{R) / (Jm{R), where (Tg{R) and am{R) are the rms fluctuations of the galaxy and 
mass density fields, smoothed with a Gaussian filter of radius R, of the galaxy 
and mass distributions, respectively. In our full analysis, we also compare our (3 
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estimates to the function f3p{k) = Cl^-^/bp{k). Here bp{k) is the bias function 
defined as bp{k) = Pg{k) /P„i{k), where Pg{k) and Pm(k) are the power spectra 
of the galaxy and mass distributions, respectively. In the case of a simple, linear 
bias model, 5g = b5m, all these definitions of b are equivalent, i.e., b^ = bp = b. 

In most cases, we find that the bias factor b = Vf-^ / (3 derived from redshift- 
space anisotropy or POTENT is similar to the large-scale value of bp{k) or b(j[R) 
defined by rms fluctuation amplitudes. Discrepancies at the 10-20% level may 
reflect non-linear dynamics as much as non-linearity of bias, since they also occur 
for unbiased models. In greater detail: (1) Measurements of (3 derived by fitting 
the Cole et al. (1995) exponential velocity distribution model to P^{k)/P^{k) 
underestimate the large-scale values of Pp by about 10-20% for all cosmological 
models and biasing prescriptions. (2) Measurements of /3 derived by fitting the 
Hatton &; Cole (1999) empirical model to P2{k) / Po{k) accurately reproduce the 
large-scale values of /3p for all cosmological models and biasing prescriptions. (3) 
Measurements of (3{R) derived from a POTENT-like comparison of — (V- v), the 
negative divergence of the peculiar velocity field, to Sg, the galaxy density field, 
underestimate the large-scale values of Pa by about 10-20% for all cosmological 
models and biasing prescriptions, with the exception of the Sqrt-Exp. biasing 
prescription. The Sqrt-Exp. model gives erratic results becaiisc the non-linearity 
of the bias model makes the estimated /3 sensitive to the range of Sg used in the 
linear fit. (4) Estimates of f3 from the anisotropy of the rcdshift-space power 
spectrum approximately agree with those from the POTENT method on large 
scales, except in the case of the Sqrt-Exp. model. 
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Figure 1. Comparison of (3 estimates from different methods, as 
a function of scale. Each panel shows f5{R) for a particular cosmol- 
ogy and biasing prescription. The dashed line represents the func- 
tion l3a{R) = ^^'^ /ba{R). Here, ba{R) is the bias function defined as 
bfj{R) = cry{R)/am{R), where (Jg{R) and <Jm{R) are the rms fluctua- 
tions of the galaxy and mass density fields, smoothed with a Gaussian 
filter of radius R. The solid and dotted lines show the estimates of /3, 
along with their la uncertainties, derived from fitting the exponential 
velocity distribution model (Cole et al. 1995) to P^{k)/P^{k). The 
points represent POTENT-like /3 estimates, derived by fitting a line 
to the measured relation between — (V • v) and 5g, when the velocity 
and density fields are smoothed with Gaussian filters of different radii, 
R. Each point represents the weighted average over four independent 
simulations, and its la error is computed as the error in the weighted 
mean. With the exception of the galaxy distributions created using 
the Sqrt-Exp. biasing prescription, the large-scale values of (3 derived 
from the — (V • v) — 6g relation are similar to the large-scale values of 
underestimating them by about 10-20%. Furthermore, these esti- 
mates are in approximate agreement with jS estimates derived from the 
anisotropy of the redshift-space power spectrum. 
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